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ABSTRACT 

We first present a self-consistent dynamical model in which ui Cen is formed from 
an ancient nucleated dwarf galaxy merging with the first generation of the Galactic 
thin disc in a retrograde manner with respect to the Galactic rotation. Our numerical 
simulations demonstrate that during merging between the Galaxy and the u> Cen's host 
dwarf with Mb — — 14 mag and its nucleus mass of 10 7 M Q , the outer stellar envelope 
of the dwarf is nearly completely stripped whereas the central nucleus can survive 
from the tidal stripping because of its compactness. The developed naked nucleus has 
a very bound retrograde orbit around the young Galactic disc, as observed for uj Cen, 
with its apocenter and pericenter distances of ^ 8 kpc and ~ 1 kpc, respectively. 
The Galactic tidal force can induce radial inflow of gas to the dwarf's center and 
consequently triggers moderately strong nuclear starbursts in a repetitive manner. 
This result implies that efficient nuclear chemical enrichment resulting from the later 
starbursts can be closely associated with the origin of the observed relatively young 
and metal-rich stars in uj Cen. Dynamical heating by the uj Cen's host can transform 
the young thin disc into the thick disc during merging. 
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1 INTRODUCTION 

The most massive Galactic globular cluster uj Cen is ob- 
served to have unique physical properties, such as a very flat- 
tened shape for a globular cluster (e.g., Meylan 1987), broad 
metallicity distribution (e.g., Freeman & Rodgers 1975; Nor- 
ris et al. 1996) , strong variations of nearly all element abun- 
dances among its stars (e.g., Norris & Da Costa 1995; Smith 
et al. 2000), kinematical difference between its metal-rich 
and metal-poor stellar populations (e.g., Norris et al. 1997), 
multiple stellar populations with different spatial distribu- 
tions (e.g., Pancino et al. 2000; Ferraro et al. 2002), star 
formation history extending over a few Gyr (Lee et al. 1999; 
Smith et al. 2000), and its very bound retrograde orbit with 
respect to the Galactic rotation (Dinescu et al. 1999). These 
unique characteristics have been considered to suggest that 
there are remarkable differences in star formation histories, 
chemical enrichment processes, and structure formation be- 
tween uj Cen and other Galactic normal globular clusters 
(e.g., Hilker & Richtler 2000, 2002) 

The observed extraordinary nature of uj Cen has at- 
tracted much attention from theoretical and numerical 
works on chemical and dynamical evolution of uj Cen (e.g., 



Icke & Alcai'no 1988; Carraro & Lia 2000; Gnedin et al. 2002; 
Zhao 2002). One of the most extensively discussed scenario 
for uj Cen formation is that uj Cen is the surviving nucleus of 
an ancient nucleated dwarf galaxy with its outer stellar enve- 
lope almost entirely removed by tidal stripping of the Galaxy 
(Zinnecker et al. 1988; Freeman 1993). The observed atyp- 
ical bimodal or multi- model metallicity distribution (e.g., 
Norris et al. 1996) and the metal-rich stellar population 2 — 
4 Gyr younger than the metal-poor (Lee et al. 1999; Hilker 
& Richtler 2000; Hughes & Wallerstein 2000) have been sug- 
gested to support this scenario. However, because of the lack 
of extensive numerical studies on dynamical evolution of nu- 
cleated dwarf galaxies interacting/merging with the Galaxy, 
it remains unclear when and how an ancient nucleated dwarf 
galaxy lose only its stellar envelope without totally destroy- 
ing its nucleus in its dynamical interaction with the Galaxy. 

The purpose of this paper is to demonstrate that uj Cen 
can be formed from an ancient nucleated dwarf galaxy inter- 
acting/merging with the young Galactic disc (~ 10 Gyr ago). 
We consider that if uj Cen is formed from merging between a 
massive, compact nucleated dwarf and the Galaxy, the merg- 
ing epoch should be well before the formation of the present- 
day thin disc, because such a massive dwarf can significantly 
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Figure 1. Orbital evolution of the nucleated dwarf for four dif- 
ferent models with initial e p ~ 0.63: The fiducial model with Mb 
= —14.0 mag, Fa =5, and Mj, = (thick solid), less luminous 
one with Mb = —13.0 mag, Fd =5, and Mb = (thin solid), 
bulge model with M B = -14.0 mag, F d =5, and M b = 10 10 M 
(thin dotted), and the more massive Galaxy model with Mb = 
— 14.0 mag, Fd =20, and M^ = (dashed). The upper and lower 
thick (horizontal) lines represent the observed apocenter (6.2 kpc) 
and pericenter distance (1.2 kpc) of ui Cen's orbit (Dinescu ct al. 
1999). Note that only the nucleus of the more luminous dwarf can 
reach the central region of the Galaxy within a few Gyr. Note 
also that the dwarfs in the latter two models can not approach 
the inner region of the Galaxy because the dwarfs are completely 
destroyed before dynamical friction cause significant orbital decay 
of the dwarf. 

heat up the thin Galactic disc (e.g., Quinn et al. 1993). Our 
fully self-consistent numerical simulations demonstrate that 
the stellar envelope of the nucleated dwarf with Mb ~ — 14 
can be nearly completely stripped by the strong tidal field of 
the first generation of the Galactic thin disc with the stellar 
mass only ~ 10 % of the mass of the present-day Galactic 
thin disc (i.e., the same as that of the present-day thick disc) 
whereas the central nucleus can remain intact owing to its 
compactness. Recently Mizutani et al. (2003) and Tsuchiya 
et al. (2003) have discussed the formation of u> Cen in terms 
of tidal disruption of a dwarf by the present-day Galaxy. 
We discuss the origin of the relatively metal-rich and young 
stellar populations of uj Cen in terms of the dwarf's star for- 
mation history strongly influenced by tidal interaction with 
the first generation of the Galactic thin disc. 



2 THE MODEL 
2.1 The Galaxy 

We construct the dynamical model of the Galaxy embedded 
in a massive dark matter halo by using Fall-Efstathiou model 
(1980) with the total mass of M t . The exponential disc of the 
Galaxy is assumed to have the radial scale length (ad) of 3.5 
kpc, the size (Rd) of 17.5 kpc, the mass of Md, and the ratio 
of the total Galactic mass to the disc mass (M t /M d , here- 
after referred to as F d ). These M d and F d are considered to 
be important parameters for the young Galactic disc. In the 
adopted Fall-Efstathiou model, the rotation curve becomes 
nearly flat at 1.75 a d and the dark matter halo is truncated 
at 8.0 aa, which is beyond the region reached by adopted 
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Figure 2. Morphological evolution of stellar components of the 
nucleated dwarf galaxy projected onto the x-z plane (upper four) 
and onto the x-y one (lower four) in the fiducial model. For clarity, 
the Galactic plane is represented as a solid line in the three of the 
upper four panels. The time T (in our units) represents the time 
that has elapsed since the simulation starts. Each frame in the 
upper (lower) four panels measures 54.6 (9.4) kpc on a side. For 
comparison, wc plot the initial Galactic disc stars at T = and 
the orbit for the first 0.93 Gyr (dotted) in the upper left panel of 
the upper four. 

orbit of uj Cen. If we assume that M d = 6.0 x 10 10 M 
(hereafter referred to as M d (0)) and M t = 3.0 x 10 11 M 
(hereafter M t (0)), the maximum circular velocity is 220 km 
s _1 . In addition to the rotational velocity associated with 
the gravitational field of disc and halo component, the ini- 
tial radial and azimuthal velocity dispersion are given to disc 
component according to the epicyclic theory with Toomre's 
parameter (Binney & Tremaine 1987) Q = 1.5. The vertical 
velocity dispersion at given radius are set to be 0.5 times 
as large as the radial velocity dispersion at that point, as is 
consistent with the observed trend of the Milky Way (e.g., 
Wielen 1977). 
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In the present study, we consider that merging between 
the lo Cen's progenitor and the young Galaxy happened 
when the stellar disc mass of the Galaxy is only 10 % of 
the mass of the present-day Galactic thin disc (M d (0)). The 
total mass of the young Galactic disc is highly likely to be 
larger than 0.1 (because of the possibly gas-rich nature of 
the disc) and we mainly present the results of the models 
with Md = 0.2 x Md(0). Since there are no observational 
constraints on the dark matter content of the young galaxy, 
we mainly investigate the models with Fd = 10.0 (M t =0.4 
x M t (0)) and 20.0 (M t = 0.8 x M t (0)). For comparison, we 
also investigate the model with a Galactic bulge with the 
R 1 ' 4 density profile, the effective radius of 0.7 kpc, and the 
mass (M b ) of 10 1() M Q . 

2.2 The lo Cen's progenitor: a nucleated dwarf 

In estimating the initial stellar mass of the nucleated dwarf 
(i.e., lo Cen's host galaxy), we consider the following two 
points: (1) The mass fraction of the nucleus (hereafter re- 
ferred to as / n ) is observed to range from 2 % to 20 % 
for nucleated dwarfs in nearby clusters (e.g., Binggeli & 
Cameron 1991) and (2) Given the observationally estimated 
small pericenter (r per i) and apocenter distance (r apo ) of the 
present-day lo Cen's orbit with respect to the Galaxy (Di- 
nescu et al. 1999), stripping of some fraction of stars through 
long-term (~ 10 Gyr) tidal interaction with the Galaxy (e.g., 
Combes et al. 1999) is highly likely for lo Cen. We therefore 
consider that the initial stellar mass m dv of its host is equal 
to (1.0 — /lost) 1 x fn 1 x m w , where is the present-day 
lo Cen's mass (= 5.0 x 10 6 M ; Meylan et al. 1995). If 
we adopt /i oat of 0.2 and / n of 0.05, m dw is 1.25x 10 s M Q . 
The nucleated dwarf with Mb (B-band absolute magnitude) 
is assumed to consist of dark matter halo, stellar envelope, 
and nucleus. The stellar nucleus is modeled by the Plum- 
mer model with the scale length of a n . Although we inves- 
tigated both dwarf disc and spheroidal/elliptical models for 
the stellar envelopes, we show only the results of the dwarf 
disc models. We use the Fall-Efstathiou model (1980) for 
the dwarf disc models with M t = 10 x Md and the central 
B— band surface brightness (no) of 22 and 24 mag arcsec~ 2 . 

2.3 The dwarf's orbit 

The center of the Galaxy is set to be (x,y,z) — (0,0,0) and 
the initial position and velocity of a dwarf are (x,y,z) = (cos# 
x n n , 0, sin# x ?"i n ) and (v x ,v y ,v z ) = (0, v ln , 0), respectively, 
where n n , 9, and, Vi n are the distance from the Galactic cen- 
ter, the inclination angle of the dwarf's orbit with respect 
to the Galactic plane, and the velocity of the dwarf, respec- 
tively. The positive sign of Vi n represents a prograde orbit 
with respect to the Galactic rotation. We present the results 
for the models with r ln — 1.5 Rd (= 26.25 kpc), in which 
the dwarf can intrude into the Galaxy from well outside the 
disc component. 

2.4 The fiducial model 

We first searched for a collisionless model in which the end 
product satisfies the following two conditions: (1) Only the 
dwarf's nucleus can survive with the envelope being almost 



completely stripped and (2) the developed "naked nucleus" 
has an orbit similar to that observed for the present-day lo 
Cen (i.e., strongly retrograde orbit nearly confined within 
the Galactic plane; e.g., Dinescu et al. 1999). We found that 
the model with Mb ~ — 14 mag, the dwarf's nucleus mass 
fraction of ~ 0.05, no = 24 mag arcsec - , 6 ~ 30°, and Vi n 
~ —60 km s _1 (i.e., e p ~ 0.63) can satisfy the above two re- 
quired conditions in the bulgeless Galaxy model with Ma = 
0.2 x M d (0) and F d = 10. We mainly show this fiducial re- 
sults of the model. Fig. 1 shows the orbital evolution for the 
fiducial model as well as for those which did not succeed in 
satisfying the two conditions. These unsuccessful collision- 
less models suggest that (1) it is the density of the satellite 
that determines at which galactiocentric radius it will dis- 
rupt, while it is the mass of the satellite that determines the 
rate of orbital decay and (2) the more massive young Galaxy 
and the Galaxy with the bulge can prevent the survival of 
lo Cen-like clusters with the pericenter of < 5 kpc. 

For the fiducial model, we also investigate the star for- 
mation history of the dwarf disc by assuming the gas mass 
fraction of 0.1 and by adopting the Schmidt law (Schmidt 
1959) with an exponent of 1.5. An isothermal equation of 
state is used for the gas with a temperature of 2.5 x 10 3 
K. We first describe the results of the collisionless fiducial 
model and then describe the model including star formation. 
The total number of particles used for each model are 30000 
for the Galactic dark matter, 20000 for the Galactic stellar 
disc, 10000 for the dwarf's dark matter, 20000 for the dwarf's 
stellar envelope, 10000 for the dwarf's gas, and 5000 for the 
dwarf's nucleus (All of these are "live"). All the calculations 
related to the evolution of collisionless models have been car- 
ried out on the GRAPE board (Sugimoto et al. 1990) and the 
models including star formation and hydrodynamical evolu- 
tion are investigated by using TREESPH codes described in 
Bekki (1995, 1997). Different gravitational softening lengths 
are allocated for different components (e.g., nucleus) so that 
we can investigate the dynamical scale on both the Galaxy- 
scale and the dwarf's scale. 



3 RESULTS 

Figs. 1 and 2 summarise the dynamical evolution of the stel- 
lar component of the dwarf in the collisionless fiducial model. 
As the dwarf sinks into the inner region of the Galactic disc 
owing to dynamical friction, the outer low surface brightness 
stellar envelope is efficiently stripped by the Galactic strong 
tidal force (T = 0.93 Gyr). The compact prolate bar-like 
structure can be formed when the dwarf passes by its orbital 
pericenter because of the tidal perturbation. The dwarf fi- 
nally lose most of its initial stellar (and dark matter) mass 
within ~ 2.6 Gyr. The stripped stars forms an inner stel- 
lar halo with the total mass of ~ 10 8 M (corresponding to 
less than 10 % of the present-day Galactic stellar halo mass) 
and the initial Galactic thin disc finally become significantly 
thickened owing to the vertical heating by the dwarf (See 
Fig. 3). 

The nucleus, on the other hand, can survive tidal de- 
struction by the Galaxy owing to its initial compact con- 
figuration (See Fig. 3). The developed "naked nucleus" still 
follows the decayed dwarf orbit at the time of dwarf's de- 
struction so that it has orbital eccentricity of 0.78 (r por i = 
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Figure 3. Upper two: Final mass distribution of the Galactic 
stellar disc (left) and the survived dwarf's nucleus (right) pro- 
jected onto the x-z plane at T = 2.56 Gyr. Note that the w 
Cen's host dwarf can cause significant dynamical heating (ver- 
tical thickening) of the first generation of the Galactic thin disc. 
Lower two: Initial (dotted) and final (solid, at T = 2.56 Gyr) ra- 
dial dependences of the vertical velocity dispersion (cr z ) for the 
initial Galactic disc in units of km s~ x (left) and the initial (dot- 
ted) and final (solid) density distributions of the dwarf projected 
onto the x — y plane for the central 200 pc. The decrease of the 
central density of the dwarf is due mostly to numerical effects 
caused by the introduction of gravitational softening length. 



1.0 kpc and r apo =8 kpc for the final 0.2 Gyr) similar to the 
observationally suggested one (~ 0.7). The total mass within 
the central 100 pc of the dwarf is 6.2 x 10 M Q for the nu- 
cleus, 3.4 x 10 6 M© for the stellar envelope, and 0.4 x 10 6 
Mq for the dark matter at T = 2.6 Gyr, which means that 
the surviving nucleus is dominated by baryonic components. 
Thus a u) Cen-like stellar cluster with the mass of 10 7 M 
and almost no dark matter can be formed from the dwarf 
galaxy dominated by dark matter. We expect that this clus- 
ter will further slowly decrease its apocenter distance and its 
stellar mass owing to dynamical friction and tidal stripping 
(e.g., Zhao 2002). 



4 STAR FORMATION HISTORY 

Fig. 4 shows that the star formation of the fiducial model 
with gas dynamics is moderately enhanced around T ~ 1.2 
Gyr (0.025 M yr" 1 ) and ~ 2.2 Gyr (0.02 Mq yr" 1 ; also at 
T = 0.2, 0.6, 0.8, 1.0, and 1.5 Gyr). This enhancement of 
star formation can result from the radial gas inflow induced 
by tidal torque of the developed prolate (bar-like) stellar 
structure in the dwarf. As a result of this, 40 % of the initial 
gas is consumed up by the triggered star formation, and 
mass fraction of new stars to the nucleus (old stars) within 
200 pc becomes rather high (0.21) at T = 2.6 Gyr. 

Since such efficient star formation can naturally cause 
rapid chemical enrichment due to metal ejection from super- 
novae, the new stars are highly likely to be more metal-rich 
than the old nucleus component and form a secondary peak 
of the metallicity distribution of the nucleus. These results 
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Figure 4. The time evolution of star formation rate (in units of 
Mq yr _1 ) in the fiducial model (left) and the time evolution of 
the mass (normalized by the initial gas mass) for the gas (solid) 
and the new stars (dotted) in the model (right). 



suggest that the observed unique chemical evolution history 
of uj Cen can be understand in the context of the repetitive 
radial gas inflow to the dwarf's nucleus (and the resultant 
moderate starbursts and chemical enrichment there) trig- 
gered by tidal interaction with the Galaxy. We also suggest 
that the observed spatial distribution and kinematics of rel- 
atively metal-rich populations with —1.2 < [Fe/H] < —0.6 
in u> Cen (Norris et al. 1997; Pancino et al. 2000; Ferraro et 
al. 2001) can reflect the detail of the radial gaseous inflow 
processes in the dwarf nucleus. The new stellar populations 
in the nucleus might well be classified as new clusters rather 
than field stars, because star clusters are more likely to be 
formed in the central regions of interacting/merging galax- 
ies because of very high gas pressure (Bekki et al. 2002). 
The merging of these young clusters with the dwarf's nu- 
cleus can cause flattening (by rotation) of the shape of the 
nucleus and thus be responsible for the observed flattened 
shape and unusual rapid rotation of w Cen. 



5 DISCUSSIONS AND CONCLUSIONS 

If uj Cen was previously the nucleus of a nucleated dwarf 
galaxy, what fossil evidences for this can be seen in the 
Galactic halo region ? Fig. 5 demonstrates that the tidally 
stripped stellar envelope of the uj Cen's host show a char- 
acteristic distribution in the L z — L xy plane, where L z (L x , 
L y ) and L xy are the angular momentum component in the z 
(x, y) direction and (L x + Ly) 1 ^ 2 , respectively. These stars 
around the solar neighborhood also show some crowding 
around L z ~ —500 and L xy ~ 300 kpc km s _1 , which reflects 
the orbital evolution of the dwarf. If we adopt the observed 
luminosity-metallicity relation [Fe/H],, = — 3.43(±0.14) — 
0.157(±0.012)M V (Cote et al. 2000) for dwarfs, we can ex- 
pect that the stellar halo formed from the u> Cen's host with 
Mb ~ — 14 mag has the likely peak value of [Fe/H] ~ —1.2 
(or somewhere between —1.5 and —0.84 in [Fe/H]) in its 
metallicity distribution for B — V = 0.5. Thus we suggest 
that the Galactic halo stars with [Fe/H] ~ —1.2 and L z ~ 
—500 and L xy ~ 300 kpc km s _1 can originate from uj Cen's 
host. 

uj Cen-like objects have been already discovered in other 
galaxies and environments: Gl in M31 (e.g., Meylan et al. 
2001) and very bright Gl-like cluster in NGC 1023 (Larsen 
2001). We suggest that if uj Cen-like objects in disc galax- 
ies are formed from ancient nucleated dwarfs merging with 
discs, there should be some correlations between the exis- 
tence of lo Cen-like objects and structural properties of discs, 
because galaxy interaction/merging can be responsible not 
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Figure 5. The distribution of stars stripped from u> Cen's host 
dwarf on the L z — L xy plane for all stars (left) and those within 
5kpc from the solar neighborhood (right). Here L z (L x , L y ) and 
L X y arc the angular momentum component in the z (x, y) direc- 
tion and (L x + Ly) 1 ' 2 , respectively, and this L xy is not strictly 
a conserved quantity. The right panel can be directly compared 
with observations shown in Fig. 15 of the paper by Chiba & Beers 
(2000). 



only for the formation of thick discs and bars and for the 
formation of starbursts and AGNs (e.g., Noguchi 1987). For 
example, it is an interesting observational question whether 
or not disc galaxies with w Cen-like objects are more likely 
to have thick discs. 
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